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Abstract 

In this paper, we derive asymptotic models for the propagation of two and three-dimensional gravity 
waves at the free surface and the interface between two layers of immiscible fluids of different densities, 
over an uneven bottom. We assume the thickness of the upper and lower fluids to be of comparable size, 
and small compared to the characteristic wavelength of the system (shallow water regimes) . Following 
a method introduced by Bona, Lannes and Saut based on the expansion of the involved Dirichlet- 
to-Neumann operators, we are able to give a rigorous justification of classical models for weakly and 
strongly nonlinear waves, as well as interesting new ones. In particular, we derive linearly well-posed 
systems in the so called Boussinesq/Boussinesq regime. Furthermore, we establish the consistency of 



the full Euler system with these models, and deduce the convergence of the solutions. 

1 Introduction 

■ 1.1 General settings 

This paper deals with weakly and strongly nonlinear internal waves in a two-fluid system. We consider the 
case of uneven bottom topography and free surface, though the rigid lid assumption is mentioned. The 
idealized system studied here consists in two layers of immiscible, homogeneous, ideal, incompressible and 
irrotationnal fluids under the only influence of gravity. 

The mathematical theory of internal waves, following the theory of free-surface water waves, has at- 
tracted lots of interest over the past decades. We let the reader refer to the survey article of Hclfrich and 
Melville [15] for a good overview of the ins and outs on this problem. 

Many models for a two-fluid system have already been derived and studied. Systems under the rigid-lid 
assumption have first been investigated (see [22] or [19] for example). Weakly nonlinear models in the free- 
surface case have been presented by Camassa and Choi [10]. Nguyen and Dias [23] presented a great deal 
of numeric simulations for such Boussinesq-type systems. Strongly nonlinear regimes have been derived by 
Camassa and Choi [9], and Barros, Gavrilyuk and Teshukov [3], generalizing the classical Green-Naghdi 
model (see [13]). A different approach has been carried out by Craig, Guyenne and Kalisch [11], using 
the Hamiltonian formulation of the Euler equations. Most of these works are formal, and restricted to 
two-dimensional flows, or to the flat-bottom case. Finally, we refer to the work of Bona, Lannes and 
Saut [7] who, following a strategy initiated in [4-6], rigorously derived a large class of models in different 
regimes, under the rigid-lid assumption. This paper is concerned with the more complex case where the 
rigid-lid assumption is removed and replaced by a free surface. 
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The strategy consists in rewriting the full system as a system of four evolution equations located 
on the surface and the interface between the two fluids (as opposed to two equations in the rigid-lid 
case). The reformulation introduces a Dirichlet-to-Neumann operator G[Q and an interface operator H[(], 
defined precisely below. The computation of asymptotic expansions of these operators leads to the models 
presented here. We focus here on shallow water regimes, allowing strongly nonlinear waves. 

Our analysis gives a rigorous derivation of most of the models existing in the literature, and also 
interesting new ones. In particular, we derive a set of models in the Boussinesq/Boussinesq regime, with 
coefficients that can be chosen so that the system is linearly well-posed. We prove that the full Euler system 
is consistent with each of our models, which roughly means that any solution of the full system solves the 
asymptotic systems up to a small error. We also prove, using energy methods together with consistency, 
that the solutions of our models converge toward the solution of the full Euler system, assuming that such 
solution exist. 

The paper is organized as follows. Section [l] is devoted to the reformulation of the full system, from 
the Euler equation to the "Zakharov formulation", written in dimensionless form. In Section [l.5[ we focus 
on the linearized system, and its dispersion relations are derived. From the asymptotic expansion of 
the operators G[C] and H[(] presented in Section [21 the asymptotic models under different regimes are 
rigorously obtained, and presented in Section 12.31 The consistency of the full Euler system with each of 
our models is proved. Then, Section [3] gives convergence results : we show that the solutions of the full 
Euler system tend to associated solutions of our models in the shallow- water limit. Finally, the links with 
different models already existing in the literature are presented in Section HJ for rigid- lid models [7] and 
layer-mean equations [9, 10]. The proof of Proposition 12.51 is given in Appendix. 




X 



Figure 1: Sketch of the domain 



Notation. We use the Cartesian coordinates (X, z), where z is the vertical variable, and X is the d- 
dimensional horizontal variable : X — x when d — 1 and X — (x, y) when d = 2. 
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The symbols V and A denote the gradient and Laplace operators in the horizontal variables, whereas 
V x,z and Ax. z are their (d + Invariable version. For p > 0, we also define the scaled version of the 
gradient and Laplace operators, namely V^- z := (y/p\7 T , d z ) T and A x z := ^A + d 2 . 

Given a surface T := {(X, z), z = ((X)}, we denote by d n the upward normal derivative at T : 

d n :=n- V x ,z, with n := — — 1==(-V<, 1) T - 

V 1 + l v Cr 

If we consider an elliptic operator P = \7 X , Z • -PVx.z, then the co- normal derivative associated to P is 

dZ-.= n-FV x ,z, 

that we simply denote d n when there is no risk of confusion. 

For any tempered distribution u, we denote by u its Fourier transform. We use the standard Fourier 
multiplier notation f(D)u, defined in terms of Fourier transforms by 

J{D)u := fu. 

The operator A = (1 — A) 1 / 2 is equivalently defined using the Fourier multiplier notation to be A = 
{l + \D\ 2 f/\ 

We denote by H s (M. d ) (or simply H s if the underlying domain is clear from the context) the L 2 -based 
Sobolev spaces. Their norm is written | • | , and simply | • | for the L 2 norm. 

Then for < T < oo, q e N, W q '°°([0, T]; H s (R d )) (or simply W q '°°H s , and L°°H S when q = 0) 
denotes the space of the functions f(t,X) defined on [0, T] x M d , whose derivative up to the order q in t 
are bounded in H s (R d ), uniformly with respect to t £ [0, T). Their norm is written | • \ w<liOOHS - 

Since it often appears, it is convenient to introduce for s and T > the space H s ([0, T]), made up of the 
quadruplets (Ci,C2,«i,«2) such that u u u 2 € W 1 > oo ([0,T];H s+7 / 2 (M d )) d , (i G W 1,oo ([0, T]; H s+3 / 2 (R d )) 
and C 2 G W Al ' oo ([0,r];i? s+5 /2(Rd)). Their norm is written | • | H ,. 

Finally, we denote by S + the planar strip R d x (0, 1), and by S~ the planar strip M. d x (—1, 0). We use 
the notation || • \\ RS for the usual norm of H S (S ± ) 7 and simply || • || 2 for the L 2 (S ± ) norm. We also for 
sgl and k G N introduce the spaces 

^ fe (5 ± ) = {/eP'(^):||/||^, fc <cx3}, 
where ||/|| ff3 , fc = Ej=o || A "~ J '^/||- 



1.2 The basic equations 

We assume that each fluid is irrotational and incompressible, so that we can introduce velocity potentials 
4>i (i = 1, 2) respectively associated to the upper and lower fluid layer. The velocity potentials satisfy 

A x , z & - in Q.\, (1) 

where £l\ denotes the domain of the fluid i at time t. Moreover, we assume the fluids to satisfy the Eulcr 
equation, and their respective density pt is constant, so that the velocity potentials satisfy the Bernoulli 
equation : 

dt^ + hvx.M 2 = ---gzmSl\, (2) 

^ Pi 

where g denotes the acceleration of gravity and P is the pressure inside the fluid. The kinematic boundary 
condition at the known bottom topography Lf, := {z = — /120 + b(X)} is given by 

d n (t> 2 = on r 6 . (3) 



3 



It is presumed that the surface and the interface are given as the graph of functions (respectively X) 
and (2(t,X)) which express the deviation from their rest position (respectively (X, hm) and (X, 0)) at 
the spatial coordinate X and at time t. The assumption that no fluid particle crosses the surface or the 
interface gives the following kinematic boundary conditions : 

dtCi = Vl + lVCiPcUi on Ti := {z = h w + &(t, X)}, (4) 

dtC,2 = y/1 + |VC 2 | 2 5„0i = y/l + \Vb\ 2 d n 4>2 onT 2 :={z = ( 2 (t,X)}. (5) 

Finally, we close the set of equations assuming that 

P is constant at the surface, and continuous at the interface. (6) 

Remark 1.1. The system is always ill-posed in the absence of surface tension (see [16]). We could add a 
surface tension term to the equations, and carry on the study as above. However, because of the smallness 
of this surface tension, the models obtained in that way would be the exact same as the ones we give on 
Section \2. 3\ Thus, for simplicity, we decide to omit the surface-tension term. 



1.3 Reduction of the equations 

In [26] , Zakharov remarked that the surface wave system can be fully deduced from the knowledge of the 
surface elevation, and the trace of the velocity potential at the surface. We extend it here for two fluids in 
the free-surface case. Indeed, if we introduce the traces 

Mi,*) fci<> + &(*,*)), and fa{t,X) := <fo(t, X, ( 2 {t, X)), 

then (f>2 is uniquely given as the solution of Laplace's equation (fTJ) in the lower fluid domain, with the 
Neumann condition Q on I";, and the Dirichlet condition 02 = V>2 on IV Then, <fii is obtained as 
the solution of Laplace's equation on the upper fluid domain, with the Neumann condition given by {5J 
d n 4>2 = 9 n (f>i on T2, and the Dirichlet condition (f>\ — on Y\. 

Following the formalism introduced by Craig and Sulem in [12], we first define the Dirichlet- Neumann 
operators : 



Gi[Ci,C2,6](^i> 2 ) := Vl + lVCil^ll^/uo+d' 
G 2 [C2,«#2 := v / l + |VC 2 | 2 9^2|,= C2 - 

We also define the following operator : 

H[Ci,C2,b}(^ui>2) :=V<t>i\ z =( 2 . 
Using the chain rule and the last definitions in the relation evaluated at the surface, we obtain 

A,/, MaflMf^M 1 ™, I" (Gl[Cl,C2,6](^ 1 ,^ 2 )+VCl-V^l) 2 Pi (r . 

^+.9(1 + 0) +2^11 2(1 + |VC 1 | 2 ) = ~7i' () 

where P\ is the constant pressure at the surface. Using again the Bernoulli equation for the upper and the 
lower fluid evaluated at the interface, we have 



il z =C 2 )+9C2 + ^\H[CuC2,b}(^2)\ 2 (1+|VC 2 | 2 ) = ~7i> 



R I J- r a. 1 IV7 / |2 (G 2 [C2,fc]^2 + VC2- V^ 2 ) 2 P2 

8^2+9(2 + - 2 wh 2(1 + | VC2 p) = -y 2 > (9) 
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where P 2 is the pressure at the interface, identical in ||5J) and ©, thanks to the continuity assumption in 
©■ 

Finally, using ([U, ©j the gradient of the equality ([7]) and a straightforward combination of §8§ and 
(HI) , we obtain the system of equations 

f StCi-Gi[Ci,6,&](^,V2) = o, 
9tC2-G 2 [C 2 ,6]V 2 = 0, 

9 t v^i +3VC1 + 2- v (IWil 2 ) - VA/i(Ci,6,&,V'i>V*) = 0, (10) 
9 t (V^ 2 - 7 ^ [Ci, C2, ^2)) + , 9 (1 - 7)VC 2 + ^V(|W 2 | 2 - j\H[a, (2, 6](Vi, V2, &)| 2 ) 

-VJV 2 (Ci,C2, Mi.W = 0, 

where 7 = — , and 

^ ^ ,/n (G 1 [Ci, C2,&](V'i,^)+VCi-V^) 2 

M(Cl,C2, b,tpi,tp2) = 



2(1 + |VCi| 2 ) 

(G 2 [C 2 , *#2 + VC 2 • W> 2 ) 2 - 7 (G 2 [C 2 , «#2 + VC 2 • [£, Ca, 



2(1 + |VC 2 | 2 ) 

This is the system of equations that we use to derive asymptotic models. 
1.4 Nondimensionalization of the equations 



In this subsection, we rewrite the system (|10p in dimensionless variables, introducing dimensionless pa- 
rameters which are crucial to study the asymptotic dynamics. We denote by a\ the typical amplitude of 
the surface deformation, and by a 2 that of the interface. A is the typical wavelength (assumed to be of 
the same order for the surface and the interface and in any horizontal direction). Finally, B is the order 
of bottom topography variation. 

We define the dimensionless variables 

^ X - * t MX) 

A hio A/v9mo B 

and the dimensionless unknowns 

Q(X) : = , tpi{X) : = 



pi 



Five independent parameters of the system are thus added to 7 



ai a 2 h w s hio a B 

7 — j £2:^7—, M:=-ry, := — , (3 := — . 
mo «io * "20 mo 



So, e\ and e 2 are the nonlinearity parameters and /i is the shallowness parameter. We also define the 
convenient notation 

a — — — — 

02 £2 

Remark 1.2. The scaling for nondimensionalization has been chosen considering the solutions of the 
linearized problem. 

We now rewrite the system in terms of dimensionless variables. First, we have to define the dimension- 
less operators, associated to the the dimensionless fluid domains : 

n a := {(X,z) G R d+1 ,e 2 ( 2 (X) < z < 1 + eiCiW}, 

n 2 := {(X, z) G K d+1 , - i + (3b(X) < z < e 2 ( 2 (X)}. 
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In the following, we always assume that the domains remain strictly connected, so there is a positive value 
h such that for all X eR d , 

l + eiCiPO-eaCaOX") > h > and i + e 2 ( 2 (X) - f3b(X) > h > 0. (11) 

o 

Definition 1.3. Let C, 2 and b G W 1 >°°(R d ), swc/i that Vt 2 satisfies pi]) , and suppose that V^ 2 G ii 1/2 (R d ). 
TTien u»t/i 2 unique solution in H 2 {Q, 2 ) of the boundary value problem 

A^ jZ 2 =0 in n 2) 

02 = -02 on T 2 := {z = £2(2}, (12) 

9„0 2 = O onT b :={z = -\+[3b}, 

we define G^ 5 [e 2 C 2 , (3b]ip 2 G H x / 2 (M d ) by 

G^\e 2 C, 2 ^b\i; 2 := -^.e 2 VC,2 • V0 2 | Z=£2C2 + <9 2 2 | z=e2C2 . 

Definition 1.4. Let now (1, ( 2 , and b G VF 1 '°°(R") 6e swca i/iai f2i and i7 2 satisfy Ul\). and suppose 
V4>i, V^2 G Lf 1 / 2 (R d ). Let 0i &e t/ie unique solution in H 2 (Q 2 ) of the boundary value problem 

in f2i, 

on T 1 := {z = 1 + eiCi}, (13) 
7T =^G^ [e2 C2^^2 onT 2 . 

Then we define G^eiCi, £2(2, /?6](V>i, V2) € Lf 1 / 2 (R d ) 6?/ 

Gi 1 ' <5 [eiCi,£2C2,/3o](V'i,'02) := -/xeiVCi • V0i| z= i +£lCl + d z <p llz=1+eiCl , 

andfr"'«[eiCi,e 2 Ci,j86](^i,^) G LT/ 2 (K d ) 

ff" ,5 [eiCi, e 2 C2,jS6](^i,^) = (V0i)| z=£2 <; 2 - 
In the following, when there is no possibility of mistake, we simply write : 

G2V2 :=G£% 2 C»,Pb]ih, 

GiW>i,V2) -G^hCi.esCa,^]^!,^), 

H{^ 2 ) := ff"' a [eiCi, £2(2, W>i, V2). 

Remark 1.5. TTie existence and uniqueness of such solutions <j) 2 and 0i are given by Proposition \2. 1{ 

Using these last definitions, it is straightforward to check that the system ([TU| becomes in dimensionless 
variables (where we omit the tildes for the sake of clarity): 

adtCi- iG 1 (Vi> 2 ) = 0, 
ftfc - jGfefc = 0, 

atVV'i+aVCi + fVdV^il 2 )-^^! =0, ; 




where 



3 t (V^ 2 - 7fl"(Vi,^)) + (1 - 7)VC 2 + f V(|VV- 2 | 2 - 7l Wi, ^)| 2 ) - /^VA/b = 0, 

(iG 1 (Vi> 2 ) + £iVCi- V^) 2 

M := 



2(1 + ^| ei VCi| 2 ) 
(7^2 + e 2 VC 2 • W> 2 ) 2 - 7(^2^2 + e 2 VC 2 • H{^ 2 )) 2 



2(l + ^|e 2 VC 2 | 2 ) 

We derive the asymptotic models from this system, corresponding to different sizes for the dimensionless 
parameters. 



G 



1.5 The linearized equation 

Linearizing the system (|14p around the rest state, we obtain 

aa t Ci-iG^[0,0,0](^,V2) = 0, 



ftC2-iG^ tf [0,0]^ = 0, (15) 
9 t V^ + aVCi - 0, 

d t (Vi^2 - 7^ M,<5 [0, 0,0] (V>i, V*)) + (! - 7)VC 2 = 0. 

Now, when the surface, the interface and the bottom are flat, we have explicit expressions for the operators 
G\, C?2 and H. Indeed, taking the horizontal Fourier transform of the Laplace equations in (JT2J) and |[Tc 
we obtain that 4>2 and <f>\ are solutions of the following ordinary differential equations : 

-fi\D\ 2 y + y" = 0. 

Then, using the boundary conditions, we deduce 



so that we have 



Then we obtain 



MX, z) = MX) C0Bh(VM|-D|«) + MX) tanh(^|£>|) smh{^I\D\z), 
G%' 6 [0, 0]^ 2 = yMD\ tanh(^)V> 2 . 



MX,z) =( — ^jj^ -tanh(^l£>l)tanh(VMpl)^2(X)) cosh(VMPk) 
+ ^(JQ tanh(^|D|) sinh(^I|D|z), 

so that we have 

G^[0,0,0](^i,^ 2 ) = ?f ^ nn (smh(VMlD|)^ + tanh(ffiv>2, 
cosh(y7i|L'|) 

and finally 

ff**[0,0,0]ty 1> V 2 ) = ^ffL -tanh(^|g|)tanh(^I|J|)V^ 2 (X). 

Using these expressions in the system (fl"5)) , we can easily calculate the dispersion relations. Indeed, the 
wave frequency u>±(k), corresponding to plane-wave solutions e*"'- x ' - * a ' ±t ) are the solutions of the quadratic 
equation 

4 |fcj tanh(01|fc|)+tanh(^|fc|) Jfc| 2 tanh(^I|fc|) tanh(^|fc|) 

W — = CJ +(1—7) = = U. (lb) 

VM 1 +7tanh(^|fc|)tanh(^|fc|) M 1 + 7tanh(0I|A;|)tanh(^|A;|) 

This equation has two strictly positive solutions (and their opposite) if and only if 7 < 1, corresponding 
to the case wherein the lower fluid is heavier than the upper one. This expression also appears in [11] 
and [24] . The figure [21 shows the evolution of the wave frequencies u>± , — ui± , as functions of the wave 
number k. We chose the parameters fi = 0.1, S = 1/3, 7 = 2/3. 



2 Asymptotic Models 

We derive asymptotic models for the system, by obtaining explicit expansions of the operators. Following 
the method of [7] , it is convenient to first reduce the problems (|12p and (fT3")) to elliptic equations on a flat 
strip. 
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Figure 2: Full system dispersion 



2.1 Flattening of the domain 

We define the mappings 



Ri 
R 2 



(X,z) 
(X,z) 



s+ 

(x,n{x,z)) 
s- 

(X,r 2 (X,z)) 



with r\(X, z) 
with r%{X, z) : 



z-e 2 ( 2 (X) 



" l + e 1 ( 1 (X)-e 2 ( 2 (X)' 

z-e 2 ( 2 (X) 
1/5- pb{X)+e 2 C 2 {Xy 



and denote their inverse 
Si 

S 2 

Introducing the (d + 1) x (d+ 1) matrices 
Pi 



S+ 

(X,z) 

s- 

(X,z) 



(X, Sl (X,z)) 

n 2 

(X,s 2 (X,z)) 



with si(X, z) := (1 + £l Ci(X) - e 2 ( 2 (X))z + e 2 ^(X), 
with s 2 {X, z) := (1/5 - /36(X) + e 2 ( 2 {X))z + e 2 ( 2 {X). 



1 



d z Si 



d z sj d d ,i 



-V x s 



1 



Mi d .i 
0i, d 1 



dzsjd 
0i. d 



■Vxs, 
1 



(17) 

where mi71 is the m x n zero matrix and the d x d identity matrix, we can transform the Laplace 
equations (|12p and (| 1 3|) into elliptic boundary value problems on flat strips. 

Proposition 2.1. Let (i, and b € W 1,co (B. d ) , such that fii and Q 2 satisfy Ul\) . and suppose Vtpi, 
V^2 € _ff 1 / 2 (R £i ). Then there exists a unique solution <j>\ G H 2 (S + ) and (f> 2 G H 2 (S~) to the following 
boundary value problems 

' V x , z • P 2 Vx,zh =0 m«S-, 

02 = V>2 on {z = 0}, (18) 

d n 4>2 = on {z = -1}, 



V x ,* • PiVx, z 0i = inS+, 

4>i = ipi on {z = 1}, 

d n 4>i = d n 4> 2 on {z = 0}, 



(19) 



s 



where d n (j) stands for the upward co-normal derivative associated to the elliptic operator involved : 

d n (f> := e d+ i ■ PV x ,z4>- 

Moreover, (f>i := (X,z) i— ► (j>i(X, Ti(X, z)) (i = 1,2) respectively solve the problems \13\) and ilty) . Thus, 
the operators G\ , G2 and H can equivalently be defined with 

G 2 1p2 = E d +l ■ x,z<h\z=Oi 

Gi(ipi,ip2) = e d+ i ■ PjV 'x,z<h\x=ii 
H(i>i,ik) = V^ii^o- 

Proof. The reduction of the problems IT3|) and IT2|) on the flat strip can be found on [17] (Proposition 
2.7). The coercivity condition is satisfied thanks to (jTTJ) and the assumptions on £1, C2 (see Proposition 
2.3 of [1]) : 

3fc>o,veeM d+1 ,e-p 4 e> -j-lel 2 . (20) 

k 1 1 

Thus, we just prove here the existence and uniqueness of the i? 2 -solutions 4>i (* = 1,2). Since for all 
g e iJ~ 1//2 (R d ), h G iJ 1 / 2 (R d ), one can easily construct a function w G H 1 (S + ) such that w\ z=1 = h and 
d n w\ z= o — g, (fT§|) and (|T5|) clearly reduce to the following problem 



Vx, z • PVx.^i = / in S 
9„0i =0 on T 2 := M d x {0} 



i=0 on Ti := R d x {1}, (21) 



where / G H- 1 (S + ) and P satisfies ([2D]). 

As a first step, we introduce the variational formulation of this problem. Let us define the functional 
space 

V := {v G iJ 1 (5+),7o(w) = on R d }, 

with 70 : H 1 (S + ) — » iJ 1/,2 (R d ) the trace operator on Ti. Since 70 is continuous, V, equipped with the 
scalar product of H 1 (S + ) and the corresponding norm, is a closed subspace of H 1 (S + ), hence a Hilbert 
space. A solution of the variational problem related to (f2Tj) is then a function u £ V such that 



Vv G V, / PVu • V-u = - / /«. (22) 

Since V = {u € D(iS + ), d = on Ti} is dense in V, a solution of the variational problem (|22[) is a weak 
solution of the problem (|2T1l . 

Now we can check that a[u, v) := J s+ PVm- is a continuous bilinear form. The coercivity of a is given 
by (p^D]) and a generalized Poincare inequality (see [2], Theorem 5.4.3). Finally, since b : v G V 1— ► — / 5+ /u 
is clearly continuous, the Lax-Milgram Theorem gives the existence and uniqueness of a solution u E V of 
and thus a weak solution of (|2"T1). Moreover, one has 



IMIifi <c , ||/L-i- 

The last step consists in proving that the solution u lives in H 2 (S + ), if we assume that / G L 2 . We 
introduce for h > 0, 

r h u(x,y,z) -u(x,y,z) u(x + h,y, z) - u(x,y, z) 
u h := (x, y, z) ^ = . 

Then u h is the solution (jSTJ) with f h = ThS ~ f and g h = ™ h = R , so that 
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Then we remark that for any v € H 1 {S + ), Vh(x, y, z) = j- J^ +h d x v(t, y, z)dt, so that 



I r 1 f k II i II I 

\\ v h\\ L 2 - J \\dxv\\ L2 dt < \\v\\ H i- 
Thus, one has thanks to the duality between H 1 and Hq, 

II . II ^ \(fh,v)\ ^ ll/ILall^llta , || hi 

< SUP ijr-r. < SUP || — < \\f\\ L 2- 

«6#o \\V\\ H i vdHl 

We finally have 

Since V is a Hilbert space, we deduce that there exists w £ V and a subsequence (uh k ) such that Uh k 
weakly converges towards w. Moreover, we know that Uh k converges towards d x u in 2?'(<S + ), so we deduce 

d x ueVc\H l . 

We prove in the same way that d y u £ H 1 , so that Axu e L 2 . Finally, thanks to (|20j) . we have 
\d 2 z u\ < \A x u\ + k\V x .z ■ PV x .zu\ = |Ajc«| + 
so that m € H 2 (S + ), and the Proposition is proved. □ 

2.2 Asymptotic expansion of the operators 

Since we are looking for shallow- water models (fi <C 1), we need to obtain an expansion of the operators 
in terms of fi. The method is the following. We first exhibit the expansion of the matrix Pi in terms of \x. 
Then we look for approximate solutions <j\ pp (i — 1,2) under the form : 

C pp = 0° + ^ + m 2 ^. 

Plugging this ansatz into (TTSJ) and (TIT)]) , and solving at each order of fi, gives the From which we can 
deduce the expansion of the operators, by computing the normal derivative of 4>i PP ■ 

Since (fTH|) is exactly the same problem as involved (in the case of the water-wave) in [1], we can directly 
apply the Proposition 3.8 to the lower fluid. 

Proposition 2.2. Let t > d/2 and s > t a + 1/2, Vijj 2 E H S+11 / 2 (R 2 ), C, 2 € H s+g / 2 (R 2 ) and b € 
H S+11 / 2 (R 2 ), such that J77]) is satisfied. Then one has 

\G 2 ip2 + ■ (haVih)\ H . < v 2 C , (23) 
|G 2 ^ 2 + /iV- {h 2 V^ 2 )~ ^ 2 V -T[h 2 ,(3b]V^ 2 \ H3 <^d, (24) 

withCj = C(^,l3\b\ HS+7/2+2j ,e 2 \C 2 \ H3+!s/2+2jl \Vip2\ H . +V3+3j ), andwhereh 2 := \~(3b+e 2 C, 2 is thethickness 
of the lower layer, and 

T[h, b]V := -^V(/i 3 V ■ V) + \ (V{h 2 Vb ■ V) - h 2 \7b\7 ■ V) + hVb\7b ■ V. 



Remark 2.3. To obtain the estimate h23\) , we use the approximate solution 

We need a higher order approximation to obtain namely <^ v ' 2 = 02 PP,:L + [i 2 <ft 2 , where <fi 2 can be 

obtained using the same method as in the following study. The Proposition^^ is then obtained following 
the path of Appendix\^ for the lower fluid (see [8] for a rigorous proof). 
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The study of the upper fluid is different from the one of the lower fluid, since we have now a nonhomo- 
geneous Neumann condition on the interface. In order to manage this, we first decompose (pi := (pi + (pi, 
where (pi is the unique solution of 

V x ,z • PiVx,z<Pi = inS+, 

<pi = fa on {z = 1}, (25) 

d n <£i = on {z = 0}, 

and (pi is the unique solution of 

Vx,« • AVx, 2 fi = inS+, 

<f>i = on {z = 1}, (26) 

d n cpi = G 2 ip2 on {z = 0}. 

Again, the system satisfied by (pi reduces to the water-wave problem (where the topography of the bottom 
would be given by e 2 C 2 ), so we introduce as in Remark 12.31 the approximate solutions 
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ft** : = - tihi(hi(^^- + (z- 1))A^ -{z- l)e 2 VC 2 • V^i), 

^app,2 ^a PP .l +fi 2^ 

It follows that Giipi the contribution on the Dirichlet- Neumann operator from (pi can be expanded as in 
the following Proposition. 

Proposition 2.4. Let t > d/2 and s > t + 1/2, V^i,C2 € i? s+11 / 2 (R 2 ), Ci € i? s+9 / 2 (]R 2 ), suc/i that 
ill]) is satisfied. Then one has 

|&V>i+MV-(/nVik)| H . <// 2 C 0j (27) 
+ mV ■ (fciV^i) - m 2 V ■ T[/n, e 2 C 2 ]VVi| H3 < /i 3 Ci, (28) 

with Cj = C(i,e2|C2| H , + 7/2+2 3 -,ei|Ci| K a+5/2+2i) (VV'i I ^+7/2+23)) ^ w/iere /ii := 1 + eiCi - e 2 C2 is 
thickness of the upper layer, and T[h,b]V is defined as in Proposition ] 2. 2\ 

The last step consists in computing the contribution on the Dirichlet-Neumann operator from <pi. We 
first define (p°^ v = <p° + fKp 1 + fJ, 2 <p 2 . It is straightforward that 

Pi=P° + ,P\ with P» := ( ^ °r ) and J* := ( ^ T ^ ) . 

where we have used the notations m ,n for the m x n zero matrix, and 7^ for the d x d identity matrix. 
Plugging these expansions into (|2l)j) . using Proposition ^. 21 and solving at each order, we get : 



At order 0(1) 



so that we have 



£0*0° = inS+, 

(p° = on {z = 1}, 

£d z <P° = on{z = 0}, 



0° = 0. (29) 
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At order 0(fi) : 

^ = on{z = l}, 

= -e d+1 ■ P^xA ~ V • (h 2 Vip2) on {z = 0}, 

which gives immediately 

1 = -h 1 V-(h 2 Vip 2 )(z-l). (30) 

At order 0(/i 2 ) : 

^92^2 _ _ !j ftiV . v ^ 2 _ 2eiV(i . _ £l ACi^ 2 ) in S+, 

<f> 2 = on {z = 1}, 

^9 Z 2 = V ■ T[h 2 , /36]V^ 2 + e 2 VC 2 ■ (ftiV^a + dVCMa) on {z = 0}, 

with the notation A 2 ■= V ■ (/i 2 Vi/> 2 ). This leads to the solution 

2 =/ ll ((/i 2 V • Wl 2 )(y - y + ~) - ^i(2eiVd ■ VA 2 + e 1 AC^ 2 )(y - i) (31) 
+ (V ■ T[h 2 ,/3b]Vip 2 + e 2 VC 2 • (hiVA 2 + exV(iA 2 ))(z - 1) 

This formal derivation of allows us to obtain the expansion of G\ip 2 , the contribution on the 
Dirichlct-Neumann operator from <j>i . Formally, we have 

« -Ma + (V ■ T[/i 2 , 06] V^ 2 - • {h\VA 2 ) - V ■ (/nerVCiA)) . (32) 



Summing this expansion with the one of Proposition ^. 41 gives immediately the expansion of the full op- 
erator G\{ipi, ij) 2 ). The following Proposition gives a rigorous statement of this fact ; its proof is postponed 
to Annex lAl 

Proposition 2.5. Let t > d/2 and s > t + 1/2, V^i, Vij) 2 € H S+11/2 (R 2 ), Ci € if s+7 / 2 (R 2 ), C2 € 
iP+ 9 / 2 (R 2 ) and 6 e # S+11 / 2 (R 2 ), sudi i/ittf flip is satisfied. Then one has 

\Gi(ih,ih)+tiAi+A2)\ H . <V 2 Co, (33) 
|Gi(^i,^i) + ti-Ai + A 2 ) - M 2 (V • Ti + V • T 2 - ±V • (/i 2 V.4 2 ) - V • (ftieiVCM 2 )) | H . < M 3 Ci, (34) 

Urtf/l Cj = C(i,^|6| ffs + 7/2 + 23 ,e 2 |C 2 | ffs + 5/2 + 2 jI e l|Cl| ff3 + 3 /2 + 2 jI I VV>1 Ijjs+7/2+23! | W 2 | JJS+7/2+2J ) > ^ 710" 

A := V • (/iiV^i), -4 2 := V • (h 2 S/ip 2 ), 
71 := T[fti, e 2 C 2 ]V^i, T 2 := T[h 2 , /56]V^ 2 . 

Remark 2.6. yls in Remark \2.S\ the proof of the estimate i33\) requires the approximate solution (f)^ pp,1 , 
with 

and the second estimate uses 



In Appendix\M (Steps 4 and 5), we give estimates on <f>\ — <fii PP , obtained thanks to the trace theorem and an 
elliptic estimate on the boundary value problem solved by <pi — § app . This leads to the desired inequalities, 
since Gi(ipi,tp 2 ) - d n <t)T P = 9„(0i ~ 
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The last expansion to obtain is the one of H^i,^), which is given by the following. 

Proposition 2.7. Let t > d/2 and s > t Q + 1/2, V^i, V^2 G F S+11/2 (R 2 ), Ci G H s+7 / 2 (M 2 ) ; C2 G 
i? s+9 / 2 (R 2 ) and 6 e i7 s+11 / 2 (R 2 ), suc/i t/iat fHP is satisfied. Then one has 



|H(Vi,^)-VVi| H , </xC , (35) 
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|%,^2) -V^ 1 - A/ V^ 1 (^ 1 +^ 2 )-|/ l 2 A 1 /; 1 -/ ll e 1 VCi-V^ij| Ha < // 2 Ci, (36) 

TOf/l Cj = C(i,/3|6| ffs + 7/2 + 2j ,e 2 |C2| ffs + 5/2 + 2 J ,ei|Cl| J/ s + 3/2 + 2j I |V^l|^ +7/2+2 3, |V^ 2 | H a + 7/2+2j)> «^ ^Wlg 

£/ie notations of Proposition \2.5\ 



Proof. The proof uses the estimates (j6"5|) and (|70|) on u := 0i — (fii 1 ^' . Indeed, we have to give an estimate 
for |Vw| z=0 | H3 , an d a trace theorem (see Metivier [20] pp. 23-27) gives for all s > 0, 



|Vu| 2=0 | ffs < Cst(||A s+1 / 2 V W || L2 + \\k s - l !*d z Vu\\) L * < ^||A S+1 /2 V ^ 



X,zHl 



Then, the estimate (j6"5]) allows to conclude : 



|Vli| z= o| ffJ < ^^,to(^ I e l|Cl| ff max{ to+ 2, S + 3/2},e2|C2| fl - m ax{ to+ 2, S + 3/2}))(^ 2 || n || /i -s + l/2 H | V | jj s+ 1 ) • 

The first estimate ([33]) follows from this relation, together with the estimates and (paS]) . 

As for the Proposition 12. 5[ the second estimate (IBU]) requires the use of the higher order approximate 
solution u := (pi — 0" pp ' 2 , and the result is obtained in the same way. □ 

Remark 2.8. Using the same approximate solution as for the expansion of Gi(^i, 1P2), we obtain an 
estimate one order less precise in fi than in i33\) and {3J$ . This loss of precision is not seen at the formal 
level and comes from the term, due to the horizontal scaling, which is necessary in order to have a 
uniformly elliptic operator. 

2.3 Asymptotic models 

The expansions of the operator we obtained allow us to derive asymptotic models from (|14j) . The frame 
of this study is limited to shallow water regimes over finite-depth fluids, that is to say : fi -C 1, and 5 ~ 1, 
but our method could be extended to deep water regimes (5 <C 1), as in [7]. As we see in Section IU we 
recover most of the models which have been introduced in the literature, as well as interesting new ones 
(the Boussinesq/Boussinesq model with coefficients and the higher order system Furthermore, 
we show rigorously that (fH]l is consistent with all of these models, in the following sense (see [6]). 

Definition 2.9. The internal-wave system ^TJ) is consistent with a system S of 2d + 2 equations, if any 
sufficiently smooth solution of {1$ such that ill]) is satisfied solves S up to a small residual called the 
precision of the asymptotic model. Throughout this paper, the precision is given in the sense of L°°H S 
norms, which means that the H s norm of the residual is uniformly bounded, with respect to t where the 
solution is defined. 

Remark 2.10. The consistency does not require the well-posedness of and only concerns the prop- 

erties of smooth solutions of the system. However, if we assume the existence of such functions, we can 
prove that they are approximated by the solutions of consistent systems, as we see in Section^ 
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2.3.1 The shallow water /shallow water regime : /i< 1 

We assume here that both layers are in the shallow- water regime (p <C 1), whereas strong nonlinearity are 
allowed (ei, e 2 = 0(1)). We use the first order expansions (|23[l . (f3"3"f and J3S|), and we plug them into (fTT 
We obtain, discarding the O(^) terms, the following system : 



' ad t d + V • (hiVipi) + V • (feWa) = 0, 

a t c 2 + v-(/ l2 vV'2) = 0, 



f2 



V(|V^| 



(37) 



+ (1 - 7)VC 2 + 7"VCi + yV (|VV- 2 | 2 ) = 0, 



where hi = 1 + eiCi - £2(2 and /i 2 = f - /?& + e 2 C2- 

Remark 2.11. This system has already been introduced in the flat bottom case in [11], and equivalently, 
though under a different form, in [10]. We say more about this in Section\ 



Proposition 2.12. The full system fi^[ ) is consistent with {Sty , at the precision [iC§, with 



c Q = ca,p\b\ 



w 1 



3ff s + 7/2 



£2 C2 



>, W>2 



W 1 



|^l,ooJJs+3/2> I Wl 1^1,00^8+7/2 1 

Proo/. Let t > d/2 and s > t + 1/2. Let £/ := (Ci, C2, Wi, W> 2 ) be a solution of l[T4]). such that (JTT]) is 
satisfied, and U € H s . It is straightforward to check that we have 



f ac\Ci + V • (ftiV^i) + V ■ (/i 2 V^ 2 ) = V • (friV^i) + V • (/i 2 V^2) 

a t c 2 + v • (/i 2 v^ 2 ) = v • (/i 2 v</> 2 ) + ±g 2 ^ 2 , 



iGl(^l,V2), 



cWi + aVCi + |V (|V^| 2 ) = peaVM, 
cW 2 + (1 - 7 )VC 2 + 7«VCi + yV (|V-0 2 | 2 ) 



(38) 



|V-0i| 2 ) + fie 2 VAf2- 



-f7V(|ff(Vi,^)| 2 



Except for dt(H(ipx, 1P2) — V^>i), the right-hand side is immediately bounded by /iCo, thanks to the 
estimates (fz?3"|) , ((33)) and (|35p . The estimate on the derivative is obtained as in the following. 
We use the study of Appendix E] : we derive (p)3"|) with respect to t on both sides and get 



M 2 V^-fth. 



a t « = 

a„(&«) = v • a t y + ^ 2 e d+1 • a t h - e d+1 ■ d t {p^ x z % 



s+, 

on {z = 1}, 
on {z = 0}, 



(39) 



We now need estimates on the right-hand side of the system. Directly from the definition of h, we have 



\dth\ 



/fs+3/2 



1 <C . 



(40) 



Thanks to the Step 4 of Section IA.2) we have 

||5 t (P^)V 



x. 



|jJs+3/2,l — Oq. 



Finally, we can obtain the estimate on dtV, using the same method as here on the lower layer : 



1 



\d t V\ HS <^C(-, 



-0- 



Then we use the study of Appendix [A"l and obtain the estimates of Steps 4 and 5 for dtu, and use them as 
in Proposition 12 . 71 in order to obtain the desired inequality : 

|9 t (#(V>i,^ 2 )-Wi)| HS = \Vd t u\ Hs < M C . 



□ 
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Conservation laws. The first two equations of ([57]) reveal the conservation of mass, since a straightfor- 
ward linear combination gives 

dthi + e 2 V • (>iVV>i) = 0, 

d t h 2 + e 2 V ■ (h 2 \7^ 2 ) = 0. 1 j 

We can play with the system to obtain other conservation laws. The conservations of total momentum 
and energy are given by 

d t {"(h\u\ + h 2 u 2 ) +Vp+ (jhi + h 2 )f3Vb + V • ( 7 /iiiti ® u\ + h 2 u 2 ® u 2 ) = 0, 



^t(^(7^i|ui| 2 + h 2 \u 2 \ 2 ) + p) + iv • (7/ii|ui| 2 ui + /i 2 |u 2 | 2 u 2 ) + V • ( 7 /i 2 Vui + /i 2 u 2 + 7 /i 1 /i 2 (ui + u a )) 

+ ( 7 /nui + h 2 u 2 )(3Vb = 0, (42) 

with the notations hi = 1 + eiCi + e 2 C2, ^2 = y — /?& + e 2 C 2 i u i = (-2^4>i (i = 1,2), and the "pressure" 
p := hlh\ + hh\ + lh x h 2 . 



Dispersion relations. When we calculate the linearized dispersion relations as in Section [TT5| we obtain 
that uj±(k) satisfy : 



•4m - Lt£±^g±^ w . 



This dispersion relation is not the same as the one of the full system (it corresponds to the first order of 
the expansion in /x of the solutions of (HU)), but we still have the condition 7 < 1, for the system to be 
linearly well-posed. The figure [3] presents shallow water/shallow water model dispersion, compared with 
the dispersion of the full system, with the parameters [i = 0.1, 5 = 1/3, 7 = 2/3. 




Figure 3: The shallow water/shallow water model dispersion 
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2.3.2 The Boussinesq/Boussinesq regime : n ~ 62 ~ t\ <C 1 

In this regime, the shallowness and the nonlinearity are supposed to be small and of the same size. This 
time, we use the second order of the expansions, and obtain 

' ad t Ci + V • (/iiVV>i) + V ■ (/i 2 V^ 2 ) = m(^AV ■ Wi + V ■ T[h 2 , f ]V^2 - ^AV • V^ 2 )), 
d t (2 + V • (ftaVV^j) = A* (V • T[/i 2 , f ]V^J , 

9 t V^ x + aVCi + I V (|Wi | 2 ) - 0, (43) 
a t V^ 2 + (1 - 7 )VC 2 + «7VCi + |V (|W> 2 | 2 ) = ^7ft(yVA^ + |VAV>i), 



with T[/i, 6]^ defined as in Proposition [ 
Remark 2.13. If the bottom is flat, then T[ft_2, xl^^ 2 * s s i m Pty gjp-VA^. 

Model with improved frequency dispersion. This model is linearly ill-posed. Fortunately, following 
[4,6], we can easily derive asymptotically equivalent models, with coefficients which can be chosen so that 
the system is well-posed. For simplicity, we assume now to be in the case of flat bottom (see [8] for the 
varying bottom case). 

We rewrite the system (f43|) with new variables : Ui := V<fii(zi) (i — 1,2). From the calculations of 
Section l2~2l we obtain 



1 ,jb 



(z) = Vi - M( V ' + - l))A^i - - !) A ^2 



<P V 2 

We then define u\ and U2 as in the following : 

ui := V^'^zi) = Wi - fxhAVtpi - /i-aiAV^ 2 , 



tt 2 := V</>2 PP,1 (22) = V^2 - fJ,-^a 2 Aip 2 , 

with zi £ (0, 1) for the upper fluid, and Z2 € (—1, 0) for the lower fluid, and the coefficients 

2 2 
01 :=zi-l G [-1,0] ; a 2 := ^ + z 2 e [-1/2,0] ; 61 := ^ + a x e [-1/2, 0]. 



We plug this into (|43[) and obtain 

' a^Ci + V • (Mi) + V • (h 2 u 2 ) + ^(^V • A« x + + ^)V • Au 2 ) = 0, 

5 t C 2 + V ■ (h 2 u2) + A*^V • Au 2 = 0, 

(l + A i6 1 A)a tUl + A i^A9 tM2 + a VCi + yV(| Ul | 2 ) =0, ( 44 ) 
(1 + M (f + i)&)d t u 2 - ^Ad tUl + (1 - 7 )VC 2 + «7VCi + y V (M 2 ) = 0. 

Remark 2.14. J/ we choose a% = — 02 = — g and &i = — i } we obtain the classical "layer-mean" model 
i61\) . introduced by Choi and Camassa in [10]. As we see below, this system is linearly ill-posed. One of 
the interests of \44\ ) * s t° °ff er a large class of equivalent models, with parameters which can be chosen so 
that the system is linearly well-posed. 

Proposition 2.15. The full system \14\ ) is consistent with \44% a t the precision \J?C\, with 

V^l,oo ffs + ll/2 7 e 2|C2|^l,o cff3 + 9/2; e l|Cl|^l,cx Jff s + 7/^7|V'(/'l| H/ i, ooffs + 11/2 , |V-0 2 | wliOoKs+11/2 ). 
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Proof. Let t > d/2 and s > t + 1/2. Let U := (Ci, C2, Vipi, V^) be a solution of ||T4"]). such that (HI]) is 
satisfied, and {/ G H s+2 . 

We first give the proof for a\ — b\ = 02 = 0, corresponding to the original system (|4"3"]) . We just have 
to plug U in (Unj), as in the proof of Proposition ^. 121 Since €2 ~ we have | jue2 VA/i | ff<i + j^VA/^ < 
^ 2 Ci. The other residuals are bounded by ^ 2 C\ thanks to the estimates (j2"4")l , ([Ml) and ([55]) with e 2 -C 1, 
and the equivalent estimates on the derivatives which arc obtained as in the proof of Proposition 12.121 

The general case is obtained when we substitute Vf/'i — fibiA'Vipi — /ij-ai AV?/>2 for u%, and V^>2 — 
11-^02^1^2 for m 2 in (jHJ). We obtain (|4"3"]) up to additional terms that are clearly bounded by fi 2 C\. □ 



Dispersion relations. As we have said previously, the coefficients can be chosen so that the system 
(HU) is linearly well-posed. Indeed, it is straightforward to check from the linearized system that u>±(k), 
corresponding to plane-wave solutions e lfc X -* w ±* ; must be the solutions of the equation 

cu 4 - A((j,\k\ 2 )\k\ V + B{{i\k\ 2 )\k\\ (45) 

with 

(1 _ p iY )(l + lS( ai +l)-a 2 Y) + 7( l _ + aa)y)(1 _ {h + l )y) + (1 _ 7)( l _ a2 y ){1 



A(Y) :=- 



B(Y) :=(1"7) 



(i-fe 1 r)(i-^y) + ^a 1 r 2 

(i- a 2 y)(l - /3xF) 



{l-b l Y)(l-^ J fiY) + q- s a l Y^ 



and the notations 

l + 2oi l + 3a 2 1 + 3&1 

«i = -2j- i a2: =^^ i /3i:= ^— 

In order to have two positive solutions of (|4"5|) , the coefficients have to satisfy 02 < —1/3, and b\ < —1/2. 
We see that the original system (|4"3"|) . as well as the classical layer-mean model (|61l) are ill-posed. However, 
there exists sets of parameters <zi, 0,2, 61 such that the generalized system is well-posed. Moreover, we 
can choose the coefficients such that the dispersions meet with the ones of the full system, at the order 
3 in /J.|fc| 2 . We present in figured] the difference between the dispersion of the full system and the one of 
the Boussinesq/Boussinesq model for three sets of parameters : a\ = b\ = a 2 = corresponding to the 
original system (|43[) . a\ — — ~, 02 = — 3 and b\ = — -| corresponding to the layer-mean system (I6ip . and 
finally a\ « 0.4714, 0,2 ~ —0.3942 and b\ = —1 corresponding to optimized parameters in (|4"4"]) . Moreover, 
we chose fi = 0.1, 5 = 1/3, and 7 = 2/3. Note that except for the last set of parameters, the system is 
linearly ill-posed, so that the computation breaks for high wave numbers. 



2.3.3 The higher order system 

We are now back in the strong linearity regime, allowing large amplitude (ei, £2 = 0(1)). But now we use 
the higher order expansions (124j) . (|34|l and (|3"6"]) . and thus obtain the strongly nonlinear model 



(46) 



d t (2+A2 = tJ,V-T 2 , 

cWi + aVCi + j V (|W>i| 2 ) = A*e a VA/i, 



CW2 + (1 - 7)V( 2 + 7"VCi + y V (| W> 2 | 2 ) = ^dtVH + 7 e 2 V(V^i • VW) 

+62VAA2 + 7e 2 VA^i 
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Figure 4: The Boussinesq/Boussinesq models dispersion error 



where we have used the following notations : 

Ai := V • A 2 := V • (JiaVV*), 

7i :=T[h u e 2 C2]^i, T 2 :=T[h 2 ,f3b}^ 2 , 
H := fci(V • (ftiWi) + V • [hiVxh) - |/iiA^i - eiVCi ■ V^i), 

■/VI . 2 ) 

j\f^ ,_ (e2VC2-V^2-V-(/tlVV'2)) 2 -7(e2VC2-VV'l-V-(felV»/'2)) 2 

Proposition 2.16. The full system JJ^p is consistent with ai </ie precision jJL 2 C\, with 



c 1 = ca,0\b\ 



H /l,DO ff ,+ ll/2 ' 



e 2 K2 



2i e l|Cl|vi/i,~,H-s+V2> I^V 7 ! 



W* 1 



= /p + ll/2 



ir 1 



11/2)- 



Proof. Let i > d/2 and s > t Q + 1/2. Let U := (Ci, C2, V^i, V^) be a solution of fH|), such that dTTJ> 
is satisfied, and {/ € H s+2 . We plug C/ in (05]), and thanks to the estimates (J31J) and (Jl|, and the 

equivalent estimates on the derivatives are obtained as in the proof of Proposition 12. 12} we can check that 
the residuals are bounded by y?C\. □ 



Dispersion relations. The linearized system is exactly the same as the one of (|4"5|) . So the system is 
linearly ill-posed, and we should derive models with parameters, to obtain well-posed systems. 



3 Convergence results 

We show here how to use the consistency results obtained in Section 12.31 to prove convergence results, 
stating that solutions of (fT4"f - if they exist - remain close to the solutions of the asymptotic model on a 
relevant time scale. In order to simplify the analysis, we focus here on the shallow water/shallow water 
model ([37]) . in the flat-bottom case ((3 = 0), so that the system can be written 

BtU + A 1 {U)d x U + A 2 (U)d y U = 0, (47) 
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with U := (hi,h2,u lx ,u ly ,U2x,U2y) = (1 + eiCi - e 2C2,| + £2(2, ea^x^i, £2<9yV'i> e 2 d x ip 2 , e 2 d y ip 2 ), and the 
matrices 



A 1 (U) := 



/ 


Ulx 





hi 











\ 







U2x 








h 2 









1 


1 


Ulx 




































7 


1 








U2x 


u 2y 




V 




















/ 



A 2 (U) := 













hi 





\ 







U 2 y 











h 2 
























1 


1 


Ulx 


Uly 




























V 


7 


1 








U2x 


U 2 y J 



We first prove that our model (|47p has solutions, since it is a symmetrizable hyperbolic system. 

Proposition 3.1. Let s > | + 1. Let C/o € i? s (]R' i ) 6 , suc/i i/iai t/iere exists h > smc/i £/iai /or a/Z X m 
R rf , £/o(^0 satisfies the assumptions 



hi>h, h 2 >h, and \u lx + u ly \ < (1 — 7)/li, | u 2 x + u 2?y I ^ ^2 ■ 



(48) 





7 


7 


JUix 










\ 




7 


1 










U 2 y 






7 M lz 





7/11 















7"iy 








jhi 















W2z 








hi 







I 





U 2y 











/>2 


/ 



Then there exists V > and a unique U € C°([0, T'); i/ s (R d )) 6 suc/i i/iai [/ satisfies [ffy and ?7(i = 0) = 
Proof. We introduce the following matrix 5, namely 



It is straightforward that S(U) and S(U)A(U, £) are self-adjoint, with A(U, () := £iAi(U) + &A 2 (U). 
Then, using the Gauss reduction algorithm, one can check that S(U) is definite positive if U satisfies (|4"5|) . 
These requirements are satisfied at time t — by Uq, and we define T as the maximum time such that 
they remain satisfied for all t < T. We know that T > thanks to a continuity argument. Then since 
we have proved that S is a symmetrizer of (|4T|) - the Theorem 7.3.3 of [21] gives T' < T such that U is 
uniquely defined on [0, T"). □ 

The last step consists in proving that the solutions of (|4T|) approximate the solutions of the full system 
(|14[) . assuming that the latter exist. We first give the the following a priori estimate, which can be found 
for example in [21] (Theorem 7.3.9). 

Lemma 3.2. Let f <E C^QO; T]; H s ), F a smooth function such that F(0) = 0. Then with 
u G C 1 ([0, T]; H s ) n C ([0, T]; H s+l ) a solution of the symmetrizable hyperbolic system 



d t U + Ai(U)d x U + A 2 (U)d y U = f + F(U), 



(49) 



which satisfies for M > 

I \W 1 -°°([0,T]xR d ) 

there are constants C(M) and K(M) such that 



< M, 



\U(t)\ Hs < Ce Kt \U(0)\ Hs +C Ce K{t ~ s) \f{s)\ Hs ds 
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Proposition 3.3. We fix 7 £ (0,1) and <5 e (0, +00). For t > d/2 and s > t + 1/2, let U £ 
C 1 ([0;T];£P) 6 H C°([0;T];H s+1 ) 6 be a solution of such that {HP is satisfied and U is bounded in 
H s ([0,T]), uniformly with respect to e\, £2 € [0,1), and /x 6 (0, ^ max ]. We denote by U := (£1, £2, Ui, U2) 
the solution of \31\ ), with the same initial values, that we assume to satisfy M. Then one has 

\U-U\ HS <fiC , 

with C Q =C{l n ,8^ max ,\u\ Usl T). 

Proof. Thanks to the consistency result (Theorem I2.12[) . we know that U satisfies (|49|) . with F = and 

\f\ HS < ^C , 

with Co = C(^) \U H »). Then, the difference between the two solutions i? M := U — U satisfies (j49)) . with 
the same / and 

F(i2") := -A^i^t/ - A 2 (R»)d y U - Ai(U)d x R» - ^(tTjS^fl". 
Taking a smaller T if necessary, one has 

\U\w^°°(lO,T]xK d ) + l^1w 1 .°°([0,T]xR<i) < 

where M is independent of ei, £2 and fi. Thus, we can apply Lemma 13.21 and one has 

<cj Ce K <- t -''>\f\ II .ds<nC(±, 1 ,5^ max ,\u\ lls ,T). 

□ 



4 Links to other models 

4.1 Rigid lid in the shallow water/shallow water case 

In [7] , Bona, Lannes and Saut presented a model for internal waves in the shallow water regime, with the 
rigid lid assumption. They showed that a nonlocal operator has to appear for d — 2 (see observations 
in [14]). This operator cannot be seen in our model (|3"7| , so that it is a purely two dimensional, rigid lid 
effect. However, we show in the following how to make it appear from (|37[) . 

Indeed, the rigid lid assumption means that t\ = 0, when £2 remains > 0, so that a = 0. The system 
(|3"T|) becomes 

' V- (fciVih)+V - (/12VV2) = 0, 
5 t C 2 + V-(/ l2 V^ 2 ) = 0, 

< 9 t V^ + |v(|V^i| 2 ) =0, (50) 

cW 2 + (1 - 7)VC 2 + ?V (|W> 2 | 2 ) = 0, 

where hi = 1 — £ 2 C 2 and h 2 = \ — fib + e 2 ( 2 . 

For simplicity, we restrict ourself to the case of a flat bottom (/? = 0), but we could do the same 
calculations with j3 > 0. We first define the shear velocity 

v := V"02 — jVipi- 

From the first line, we deduce : 

V • (/law) = -V • ((hi + = ~^^V • ((1 + — - x 8e 2 C 2 )Vi>i). 

7 + 

Then we define the nonlocal operator : 
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Definition 4.1. Assuming that £ € i°°(R <1 ), rae define the mapping 

U ^ H-C 1/ 

where V is the unique gradient vector in L 2 (R d ) d , solution of the equation 

v ■ ((i + ov) = v • w. 

So from the definition, we have 

V^i - Q[^<5e 2 C 2 ]( 8 —Mv). 

7+0 7+o 

We plug this expression into (j50|) , and obtain immediately 

dtb + ^ V ■ (/ ll O[^ ( 5e 2 C 2 ] M) - 0, (51) 
ftt, + (1 - 7 )VC 2 + ?V (|« - ^jQ[^oe 2 C 2 ](ft 2 «)| 2 - ^IQ^fc^KM)!^ = 0, 



where /ii = 1 — e 2 £ 2 and hi = h + e 2 £ 2 - This is exactly the system derived in [7]. 

Using the same method, we could derive rigid-lid models from (|44p and (|46[) . The rigid- lid model in 
the Boussinesq regime has already been exhibited in [7] , and a fully nonlinear model is presented in [9] . 

4.2 The layer-mean equations 

In the literature, the water-wave system is often given by layer-mean equations (see for example [25]), 
using as unknowns the depth-mean velocity across the layers : 



1 



1+eiCi 



Ui(X):=— I V<t>i(X,z)dz with /i! := 1 + £l Ci - e 2 C 2 , 



C2C2 
E2C2 



^ /-t 2 ^2 ^ 

u 2 {X) ■=— V0 2 (X, z)dz with h 2 := - - (3b + e 2 C 2 , 

"2 J-l/8+Bb 



The systems under this form (obtained for example in [10] and [3]) are equivalent to the the system we 
derived, since one can approximate u\ and u 2 thanks to our previous unknown and tp 2 (as we see in 
the following Proposition), and conversely. Thus, our study gives a rigorous justification of these models, 
and we are able to offer consistency results. 

Proposition 4.2. Let t > d/2 and s > t + 1/2, Vt/>i, Wtp 2 G H S+11 / 2 (R 2 ), Ci G H s+7 / 2 (R 2 ), ( 2 € 
H s + 9 / 2 (R 2 ) and b G £F +11 / 2 (R 2 ) } such that (Q7J) is satisfied. Then one has 

|ur- W> x | ff , +1 <mCo, (52) 

|u 2 - V^Ijjs+i < M^o, (53) 

|Bi - VVi - /iPi(V^i,V^ 2 )| ffa+1 < v 2 C u (54) 

|u 2 - VV-2 - fiV 2 Vi; 2 \ Hs+1 < /i 2 C*i, (55) 
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with Cj = C(i,^|6| ffs+7/2+23 -,e2|C2| ffs+ 5/2 + 2 JI e i|Ci| ff3+3 /2 + 2 J - »|V^i|_ Hl>+7/a+2j , |V?/> 2 | Hs + 7/2 + 2i ), an d where 
T>i and T> 2 are defined by 

Bi(V^i,V^) = "^( r i - \{h\VA 2 ) - (/iieiVCi^)), 

A2 



with the notations of Proposition \2.5\ 

Proof. Using the Green formula with <pi the solution of (fT5]) . and a test function (p := (X, z) i— > we 
have 

-■ V- ,/ r,/- I yi' , ,. . 



<pA^- z <pidXdz = — / V^- z 0i • z (pdXdz + / ipd nx (j)\drii + / fd^fiidn 



Thus, we deduce 



Identically, we have 



= -/i / / VfadzdX + / cpG 1 (ipi,ip 2 )dX - / (pG 2 ip 2 dX. 

V-(fei!ti) = — (Gi(V'i,V'2)-G 2 V2). (56) 

V • (/i 2 u 2 ) = — G 2 V> 2 . (57) 
A* 

We now prove the estimate (|52j) . and the other are obtained in the same way. 

Using the Propositions 12.21 and 12.51 together with (|5o| . and since ([TO is satisfied, one has immediately 

|V-(Si-V^i)| HJ <fiC a , 

so that we only have to obtain an L 2 -estimate on u\ — VV>i- Using the definition of ui and the mappings 
defined on Section O we obtain 



Jo 



with^i : (X,z) e5 + n 4>i(X, si(X, z)). We deduce |ui-V^i] 2 < C(ei|Ci| W i,»o> e 2 1 Ca | ^-i 



'lift- 



■0i)||. The estimate follows now from Step 3 of Section lA.2[ together with the estimates and (|65|) . □ 
4.2.1 The shallow water /shallow water regime : fi <C 1, e = O(l) 

We use (|52|) and (f53|) in the system (|37|) . and with a straightforward linear combination, we obtain 



' dthi + e 2 V • = 0, 

d t /i 2 + e 2 V ■ (h 2 u 2 ) = 0, 

cWi + V/i 2 + /3V6 + V/ii + (|ui| 2 ) = 0, (58) 



2 

<9 f u 2 + V/i 2 + (3Vb + 7 V/n + y V (|u 2 1 2 ) = 0. 
Proposition 4.3. The full system \1J$ is consistent with f58\) . at the precision /zCo, with 

CO = C(h, P\b\ W l,oa H s + 7/2,£2\(2\ wliaoHs + 5/2 ,ei\(l 1^1,00^8+3/2 1 | V'01 | ^l.oofl-s+7/2 J | V?/>2 J ^1,00^3 + 7/2)- 

Proof. We know that from Proposition 12 . 121 that (|14p is consistent with |3T|) . at the precision /iCty. From 
and ([53"]) . we deduce that (Ci, C2, Wi, u 2 ) satisfies ([55)) up to a residual of the same order. □ 
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Remark 4.4. Note that the first two equations of \58\) are equalities (where the last two equations are first 
order approximations in fi), as we can see from i56\) and J57| j. They reveal the conservation of mass. 

Conservation of momentum and energy are the one obtained in Section \2.3.1[ when we substitute Hi for 
Vipi (j — 1,2). These conservation laws, and the one of higher order systems, had already been introduced 
in the flat-bottom case in [3]. 



4.2.2 The Boussinesq/Boussinesq regime : // ~ e 2 ~ ei <C 1 

We now restrict ourself to the flat-bottom case, since it considerably simplifies the notations, but the 
following could be derived with 0^0 without any difficulty. The estimates (l54"|) and |55|) with £2 ~ M an d 
(3 = give the following formal relations 



6 2d 

U 2 « V^2 + /U^VA^2- 



Plugging this into (|4"3")) we obtain the system 



(59) 
(60) 



' dthi + e 2 V • (h r ui) = 0, 
d t h 2 + 62V • (h 2 u 2 ) = 0, 
t ui+aVCi + yV (l^l 2 

9 t u 2 + (1 - 7)VC 2 + «7VCi + yV (|u 2 | 2 



(61) 



Remark 4.5. T/iis set of equations had been revealed in [10]. It corresponds to j44\) > with the choice of 
parameters : a,i = —\, = —\, hi = —\. This particular choice of parameters leads to a linearly ill- 
posed system. That is why it is interesting to obtain, as in Section [2. 3.21 a larger class of models, allowing 
well-posed systems. 

Since this system is a particular case of the Boussinesq/Boussinesq model (|44[) . we can apply the 
Proposition 12.151 



Proposition 4.6. The full system fi^| ) is consistent with Wl\) . at the precision fi 2 C\, with 



Ci = C(j-,[3\b\ 



e 2 <2 



W l,oo ffs + ll/2 



2 1 ^1,00 #3+11/2 



4.2.3 The higher order system 

We now do the same study, without assuming any smallness on ei, e 2 . We plug (|54"|) and ((55)) into (|4"6"|) , 
and obtain 



<9 t /l2 

<9 t ui 

9*1*2 



e 2 V • (hiui) = 
e 2 V • (h 2 u 2 ) = 



fv(N 2 ) 



/xe 2 VM + /^£ 2 V(mi • X>i) + (id t T>i, 
(1 - 7)V( 2 + 7"VCi + y V (|TT 2 | 2 ) = ufafrVH + 2?a) + e 2 V{yui • Vft 

+U 2 -X>2 +JV2+7M)), 



with the notations of Proposition 14.21 and |4"B1 when we substitute Ui for VV>i (i = 1,2). 



Proposition 4.7. The full system \1J$ is consistent with 



Ci =C(i,/3|6| 



£ 2 K2 



jyl,oo£fs + 9/2 ) 



eilCll^l.oo^s + 7/2, 



at i7ie precision [i 2 C\, with 



W l,oo ffs + ll/2 



, V^2 



(62) 



14/1,^^-5+11/2; 
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Proof. Let to > d/2 and s > to + 1/2. Let (£1, C2, V^i, V^) be a sufficiently smooth solution of (dU, such 
that (fTTj) is satisfied. We know from Proposition ^. 161 that (£1, £2, VV'i, V^a) satisfies P5|) up to a residual 
bounded by ^ 2 C\. Then, the estimates (j54"|) and |55|) give that (Ci, C2, "i, u 2 ) satisfies (jl)2"|) up to a residual 
of the same order. □ 

A Proof of Proposition 12.51 

Our proof contains three parts. First we introduce u the correction to the expansion of 4>i formally obtained 
in Section l2~2"j and we present the system solved by u. Then, we use the elliptic form of the operator to 
obtain H s estimates on u. Finally, we use these estimates to prove the desired inequalities. 

A.l System solved by u 

We first define the second order correction to the formal expansion : 

u := fa - V>i + n hi (hi l) + (z- 1)) Afa - (z - l)e 2 VC 2 ■ V^i + (z - 1)V • (/i a V^)) ■ 



From the computation carried out in Section [2.21 we know that u satisfies the following equalities : 

Vx, z ^iVx,zii = /i 2 V'PV inS+, 
u = on {z — 1}, 

d n u = G 2 02 + MV • {h 2 ViP 2 ) + M 2 «V ) on {z = 0}. 

Moreover, we notice that (|57|) gives G 2 ^ 2 + • (/i 2 VV' 2 ) = V • V, with V = ^/i 2 (VV' 2 —U2). Thus, using 
the definition of Pi in (fTT|) , we finally have the system 

V^,-P"V^« = /i 3 V^-h inS+, 

u = on {z = 1}, (63) 

d n u — V • V + n 2 ed+i h on {z = 0}, 

where we have introduced the notation V^- := (y / 7iV' r , <9 z ) t , and with h := P^V^ z 1 and 

We now give the useful estimates of the right-hand side of the system. It is straightforward to check 
that 

|p|| Hs+ i /2 ,i < C( — ,ei\(i\ Hs+3/2 ,e2\(2\ Hs+ s/2,f3\b\ Hs+s / 2 , | v"i/>i | ffs+7/2 , | VVj 2 | Hs+7/2 ). (64) 
Moreover, the Proposition ^. 21 gives 

\V\ HB < Ai 2 C(i,/3|6| ffs+5/2 ,e 2 |C2| ffs+3 /2, 1 Vi/>2 1^+5/2)- (65) 
A. 2 i7 sl -estimate (s > 0) on it 

We follow the sketch of the proof of Proposition 3 in [7], which contains five steps. 
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Step 1. Coercivity of the operator. Since £i, £2 6 W 1,oa and satisfy (fTTj) . we can check (see Proposition 
2.3 of [1]) that for any 6 € M d+1 , 

e-p^e > y\e\ 2 , 

k 1 1 

with k = ll^ijl^ + 7^(1 + m||Vsi || )• The operator is uniformly coercive in fi. 

Step 2. Existence and uniqueness of the solution. The result is given by the coercivity of the operator. 
From the assumptions on Ci, C2, b, ipi and V2, we know that h G H 8 ^ 2 ' 1 (S + ) d+1 and V G H s+1 (R d ). 
For s > 1/2, the proof of Proposition 12.11 works for the system ([6"3"|) . so that we know that there exists a 
unique solution in H 2 (S + ). We know prove by induction that for k G N, 

h G if fc+1 and V G H k => u G H k+2 . (66) 

We assume that h G H k + 2 and V G H k+1 . We thus know that u G if fe+2 , so that v := Au G if fe+1 C if 1 . 
Hence, w is the classical solution of 

V£, z -f^V^ = M 2 V^-h inS+, 

u = on{z = l}, (67) 

<9„u = V • AV" + fi 2 e d+1 ■ d x h on {z = 0}, 

with /i 2 h = /i 2 Ah + [A, P^]V£ u. Thanks to Theorem 6 of [18] : for i > f , one has 

II [A, ^]V^u|| 2 < C to ||VP"|| ||V^u|| a) 

so that h G H k+1 and Ay G if fe . The inductive hypothesis are satisfied, so that we know that v G H k+2 . 
Finally, we use the coercivity of the operator (Step 1) with the nth derivative of (p33")h and obtain 

\\d 2 d n u\\ 2 < k\\v x<z ■ P^ Xz d n u\\ 2 + \\A x d n u\\ 2 . 

It follows that u G if +3 , and is proved. The interpolation theory leads to the final result : for 
s > 1/2, there exists a unique solution u G P s + 3 / 2 of (f63l) . 

Step 3. L 2 -estimate on V^- 2 u. We multiply ()63[) by u, integrate by parts on both sides, and use the 
boundary conditions to finally obtain 



/ V^ti • P^ x . z u = [i 2 f V x<z u ■ h + f V 

JS JS J{z=0} 



u ■ V. 



From the coercivity and the Cauchy-Schwarz inequality, we deduce 

H V £, Z M ||2 ^ fe (^ 2 H h ll2ll V X,2 U ll2 + \ V \m/*\ Vu \z=o\H-m)- 

Then, a trace theorem (see Metivier [20] pp. 23-27) gives 

|Vu|^o| K _ 1/a < Cst(||Vu|| a + ||A- 1 «»Vu|| a ) 
<Cst(-^ + l)||V^u|| a . 

This finally gives the estimate 



|V^,«|| 3 <C(i ^\Ci\ wl ,^,e2\C2\ wl ,J(^H 2 + ^^\V\ H1/2 ). (68) 
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Step 4. L 2 -estimate on A S V^- z u (s > 0). We define v = A s u. Multiplying ((Ml) by A s on both sides, 
one obtains 

^x.z ■ P ^x.z v = M 2 ^x.z ■ h in K rf x (0, 1), 

i> = on {2 = 1}, (69) 

d n v = V • A s V + fi 2 e d+1 • h on {z = 0}, 

with /x 2 h = yLt 2 A s h + [A s , P^JV^ Z u. We can use Step 3 with v and obtain 

ll V *,,"][ 2 ^ ^(^^i|Ci| W r 1 .oo,e a jC| wrl . HO )( A * 2 ||A-h|| 2 + ||[A-,i=»"]V^,«|| a + ii^l^l^.^). 
We obtain the commutator estimate thanks to Theorem 6 of [18] : for s > — | and £q > f j one has 

||[A S ,/] 5 || 2 < C a , ta ||V/|| lrB1 „ {ta .._ I} ||fl|| H ._ 1 . 

In our case, it gives 

||[A 5 ,PlV^u|| 2 <C s , i0 (i e 1 |Ci| HfflK{t0+2 , s+1} ,e 2 |C2| i , niax{t0+2 , s+1} )||A s - 1 V^|| 2 . 



We finally get an estimate on ||A s V^- 2 m|| 2 in terms of ||A S 1 V^- z m|| 2 . Step 3 is the case when s = 0. By 
induction, and interpolation when s € (0, 1), we obtain the following relation for all s > : 

||A S V^|| 2 < C s , t0 (~ e 1 |Ci| ffmax{to+2 , s+1} ,e 2 |C2| Hmax oo + - + i})(MlA s h|| 2 + ^|^|y| Hs+1/2 ). (70) 
Step 5. L 2 -estimate (s > 0) on A s 0#V^ u. The equation (foil)) gives the formula 

I 1 2 | 1 2 

+ 7 51 1 d 2 z u = /i 2 V^ • h - /iV • (/iiVu - Vsi<9 z u) + nd z {Vs x ■ Vit) - d z ( + ^ 51 ! )d z u, 
from which we deduce 

||A S 9 Z 2 U || 2 < C(± ei|Ci| H , 1 ,oo,e 2 |C2| H , 1 , 00 )(/i 2 ||A s V^ z • h|| 2 + ^|| A s+1 V^ z u|| 2 ). 
Thus, we have the estimate 

||A^V^,«|| 3 <C(i, e i|CiL 1 ,., e2 |C2U 1 ^)(M 2 ||A s V^ z -h|| 2 + ^11^+^^^), 
and Step 4 allows us to conclude 

||A S <9 Z V^ Z W|| 2 < C Sito (-,ei|Cl| ffm ax{t + 2,s + 2},e2|C2| ff max{t 0+ 2, S + 2},Ai maX )(M 2 ||h|| ii - s + i, 1 + | ^ | ) • 

A. 3 Proof of the inequalities 

To obtain the first estimate, we remark that 

GiWm.V*) + + = d n u lz=1 - y 2 | ei VCi| 2 (/iiA^i - e 2 C 2 ■ VV>i + V • (/i 2 V^> 2 ) 
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We have immediately 

\ u o\ H s < C (^P\ b \ H s + l, e 2\(2\ H s + li e l\(l\ H s + l, |V^l| Ha+2) |V-02| ffs + 2 ), 

so that we just have to bound |9 ra U| z=1 |„ s . We now use the trace theorem to get 

\d n u lz=1 \ H3 < C(-,ei|Ci| w , 1 , oo! e 2 |C2| w ,i, 00 )(M|V'U| z=1 | J?8 + \d z u {z=1 \ HS ) 

< C(^,e 1 \( 1 \ wl ^,e 2 \C 2 \ wl , ac )(^/Jl\\^x,z u \\H^/^ + ll^ V x^ M ll^-i/2,o)- (71) 

The estimates obtained in Steps 4 and 5, together with (|64[) and (f6"S")) . give immediately the desired result. 

To obtain the second estimate, one has to carry on the proof with the higher order approximate solution 
obtained in Section 12.21 : 

u := (j> x - c/)^ pp ' 2 , 

and one would obtain the estimates exactly as above. We omit this technical step. 
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